The procccdings of the Physics of Accreting Compact Binarios 



1 



Symbiotic Novae 



Joanna Mikolajewska 

N. Copernicus Astronomical Center, Bartycka 18, 00-716 Warsaw, Poland 



Abstract 

The symbiotic novae are thermonuclear novae in symbiotic binary systems - 
interacting binaries with evolved red giant donors, and the longest orbital pcri- 
ods. This paper aims at presenting physical characteristics of these objects and 
discussing their place among the whole family of symbiotic stars. 

1. Introduction 

Symbiotic stars are interacting binaries involving an evolved giant (either a 
normal red giant in S-typcs or a Mira variable surrounded by an optically thick 
dust shell in D-types) transfcrring mass to a hot and compact companion, usually 
white dwarf. A typical symbiotic binary is embedded in a circumstellar nébula 
which is mainly formed from material lost in the red giant, while the hot com- 
poncnt is responsible for its ionization. The nature of the giant determines the 
orbital scparation at which the symbiotic interaction occurs: the binary must 
have enough room for the red giant (and in the case of D-types also for its dust 
shell) and yet allow it to transfer sufñcient mass to its companion. As a result, 
orbital pcriods for the S-types are of about 1-15 years, and more than 20 years 
for the D-types, which are the longest orbital periods among interacting binaries. 
Symbiotic stars are important tracers of late phases of stellar evolution, and a 
very attractive laboratory to study various aspeets of interactions and evolution 
in binary systems [3] . The presence of both the accreting white dwarf and the red 
giant (with its degencrated core) makes also symbiotic stars a promising "factory" 
of Type la supernovae (SN la), whatever the path a thermonuclear explosión of 
CO white dwarf upon crossing the Chandrasekhar limit (the single degenerate 
(SD) scenario) or by merging of a double CO white dwarf system (the double 
degenerate (DD) scenario) may be. 

Symbiotic novae form a small subclass of thermonuclear (TNR) novae which 
oceur in symbiotic binary systems. The aims of this paper are to present the 
state-of-the-art in understanding of these objects, and to discuss their place among 
symbiotic stars. 

2. The Nature of Symbiotic Hot Components and their Activity 

The hot components of symbiotic binaries were discussed by Mikolajewska ¡16} 
Q21 El]- Based on their activity, all symbiotic stars can be split into two subclasses: 
ordinary or classical symbiotic stars (majority) and symbiotic novae. 

The typical (quiescent) hot components of classical symbiotics appear to be 
quite hot (~ 10 5 K) and luminous (~ 100-10 000 L Q ), and they overlap in the 
same región in the HR diagram as central stars of planetary nebulae ([MJ ITS] , 
and Fig ÍL]) . In most cases, the luminosity is too high to be powered solely by 
accretion. However, if the symbiotic white dwarfs burn hydrogen-rich material as 
they acérete it, an accretion rate of the order of lO _8 M0/yr is then sufñcient to 
power the typical symbiotic component with Mwd ~ 0.5 Mq ¡16j . 

Many of the classical symbiotic systems, including the prototype Z And, CI Cyg 
and AX Per, show occasionaJ 1-3 mag optical/UV outbursts on timescales from 
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Fig. 1. The hot components in the HR diagram. The luminosities, and effective tem- 
peratures are derived from optical and UV emission line fluxes (Hell and Balmer 
HI) and UV continuum, assuming case B recombination and blackbody energy dis- 
tribution. Filled and open circles represent the S and D-type systems, respectively 
(Mikolajewska, in preparation). Diamonds correspond to measurements by Mürset 
et al. [24] whereas stars are data from |17] , and references therein. Lines connect 
the points corresponding to Z And-like systems during quiescence and outbursts, 
respectively. The solid curves with dots are steady models of Nomoto et al. |25j . 



months to years, when the hot component luminosity rcmains roughly constant 
whereas its effective tempcrature varies from ~ 10 5 to <~ 10 4 K (FigITTJ). The 
mechanism of this activity was a mystery for a long time, and only recently it has 
been q ualitati vely explained by unstable disc-accretion onto H-shell burning white 
dwarf [16l[3Q]. Unfortunately, it is still waiting for a good quantitative model. 

The thcory prcdicts that the luminosity of TNR-powered white dwarf should 
be rclated to its mass. Unfortunately, there is no unique core mass-luminosity 
(ML) relation for accreting white dwarfs because it dcpcnds significantly on their 
thermal history: generally, hot white dwarfs should have larger masse s th an the 
cooler ones to reach the same luminosity during the H-burning phase (¡llj). The 
positions of symbiotic white dwarfs in the luminosity versus mass plañe together 
with different ML relations are shown in Fig (5] The symbiotic white dwarfs cluster 
around the ML relations for stars leaving the AGB with the CO core and the RGB 
with the He core for the first time, and this may indicate that they could be still hot 
at the onset of the mass transfer from the red giant. The luminosities of the white 
dwarfs with masses > 0.5 M Q are generally consistent with those predicted by the 
steady models of Nomoto et al. ¡25] , however these models predict systematically 
highcr effective tcmpcraturcs than those derived from observed UV continuum 
and UV/optical emission lines (see Fig JTTj) . 

The symbiotic novae are thermonuclear novae. These are either extremely slow 
novae, their outbursts can go on for decades, or very fast recurrent novae with very 
short timcscalcs of their outburst, ~ scveral days, and recurrence time of ordcr 
of several years (Figj47J). The term "symbiotic nova" (SyNe) is usually reserved 
for the former whereas the latter are called "symbiotic recurrent novae" (SyRNe). 
The diffcrcnccs in outburst behavior seem to reflect different white dwarf masses: 
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Fig. 2. Mass-luminosity relation for the 
symbiotic white dwarfs. Systems with 
Z And-type activity are represented by 
filled circles with arrows denoting their 
luminosity changes. Non-eruptive Sys- 
tems and SyNe during the plateau phase 
are plotted as open circles, and squares, 
respectively. Accretion-powered sys- 
tems (stars) including the SyRNe at 
quiescence are also shown. The ML re- 
lation for accreting cold white dwarfs, 
the Paczyriski-Uus relation for ABG 
stars, and that for He white dwarfs QT| 
are plotted as solid, dashed, and dotted 
curves, respectively. The shaded región 
represents the steady models from [25] . 



very high in the SyRNe, and much lower in the SyNe. 

Such novae are very rare among symbiotic stars. Although over 200 symbiotic 
stars are known [2J, there are only ninc SyNe occurring in both S-type systems 
(AG Peg, RT Ser, V1329 Cyg, and PU Vul) and D-type systems (RR Tel, V2110 
Oph, V1016 Cyg, HM Sge, and RX Pup), and only four SyRNe (RS Oph, T CrB, 
V3890 Sgr, and V745 Seo), all of them related to the S-types. The most recent 
symbiotic nova V407 Cyg (2010) shares characteristics of both groups: it has a 
Mira companion but its very fast outburst dcvelopmcnt is typical for the SyRNe. 

3. Binary Parameters and Mass Transfer Mode 

The distributions of the orbital parameters and the component masses have 
been recently discussed in [TH1 HH]. Since then, new orbital periods [5ül |ü] and 
component masses [101 [7] have been derived. The updated distributions of the 
orbital parameters are shown in FigjXJ and they do not affcct the previous con- 
clusions. In particular, all but one (R Aqr) systems with measured orbital periods 
belong to the S-type, and most of them have P or b ~ 300-1000 days. 

At present, the orbital periods are known for all SyRNe, and the four SyNe 
in S-type systems. The periods of all SyRNe are below ~ 600 days whereas the 
SyNe have P or b > 800 days, consistent with the shorter and longer period tail 
of S-type systems, respectively. In particular, RT Ser and PU Vul with P or b of 
4500 and 4900 days, respectively, are among the systems with the longest known 
periods (the third one is an accretion-powered CH Cyg). It is also remarkablc 
that whereas only ~ 20% of known symbiotic binaries belong to the D-type, six 
of ten known SyNe have a Mira companion. This may be related to lower mass 
transfer and accretion rate in the SyNe than in most of classical S-type symbiotic 
stars, and the SyRNe. 

The masses of the symbiotic red giants peak around 1.6 M©, and most white 
dwarfs have masses between 0.4-0.6 M©. There are, however, significant differ- 
ences between the SyRNe and the other symbiotic binarios. In both, RS Oph and 
T CrB, the giant is the less massive component, with mass of only 0.4 — 0.8 M©, 
whereas their white dwarfs, ~ 1.1— 1.4 M©, are the most massive, sufficient to 
become SN la. The red giant envelope mass in both SyRNe is very low, less 
than ~ 0.35 M© in RS Oph, and below ~ 0.2 (and maybe cven 0.1) M© in T CrB. 
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Fig. 3. Orbital param- 
eters of symbiotic 
stars. Shaded regions 
denote locations of 
the SyNe (lighter) 
and the SyRNe 
(darker) , respectively. 



Stripping this envelope by the supernova ejecta would leave a low mass, < 0.5 Mq, 
hclium white dwarf only. So, the SyRNe with their giant donors can also account 
for the single low mass white dwarfs. 

It is usually assumed that the symbiotic binary components do not fill their 
Roche lobes, and interact via stellar wind. However, during the last decade, 
cllipsoidal variability has been detected in near infrared and red light curves of 
over a dozen classical symbiotics suggcsting that the giants are at least very cióse to 
filling their Roche lobes [TB1 [TÜ] . Such an ellipsoidal modulation is also present in 
all SyRNe but RS Oph ¡181 126] . Unfortunately, to see this modulation in classical 
symbiotic systems we need red/near-IR photometry which is hardly available. 
Thus it is possiblc that tidally distortcd donors and Roche-lobc overflow are quite 
common is symbiotic binaries with P rb below 1000 days, and the mass transfer 
and accretion is efficient enough to power the hot components of the classical 
symbiotic systems via stable/quasi-stable H-shcll burning. 

In the case of the SyRNe, relatively high, of order of 10 _7 MQ/yr, mass accre- 
tion ratcs are required to account for their short recurrence time, and the hot 
component activity between the TNR nova outbursts [TJ [8], which is easier to 
achieve by Roche-lobe overflow than by wind interaction. On the other hand, 
the slow SyNe should acérete at lower rates, below ~ 10~ 8 M Q /yr, as in the case 
of wind accretion. This explains the higher occurrences of the SyNe among D- 
type systems, and S-type systems with the longcst orbital pcriods. Moreover, the 
SyNe AG Peg is the only symbiotic binary with very good near- IR light curve and 
-forb < 1000 days which does not show any evidence for ellipsoidal variability [TÜ] , 

4. Outburst Evolution of Symbiotic Novae 

FigHDprescnts light curves from the ongoing outburst of the symbiotic nova PU 
Vul (1979), and of the two recent outbursts of the symbiotic recurrcnt nova RS 
Oph, respectively. Both light curves cover the same interval of time to show dra- 
matic differences in the outburst development between the two kinds of symbiotic 
novae. 

The very fast outburst evolution of RS Oph, and other SyRNe is due to the 
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Fig. 4. Optical and UV light curves of the symbiotic nova PU Vul (left) and the 
AAVSO optical light curve of the symbiotic recurrent nova RS Oph (right) covering 
identical time intervals. For PU Vul, the best fit model with 0.6 Mq white dwarf 
[13] is also shown. 



presence of the very massive white dwarf (e.g. ~ 1.3 M Q inRSOph[3]) accreting at 
a rate of order of lO _7 M0/yr. The luminosity of the nova at máximum is super- 
Eddington, and the constant-luminosity phase (plateau) lasts ~ 50 days only. 
The outburst material is ejected with a high initial velocity, i; exp > 4000 km/s. 
Relatively low mass of the nova ejecta, ~ 1O -7 M0 in RS Oph [31], comparable 
to the average accretion rate between the TNR nova outbursts, suggests that the 
white dwarf mass can grow. The outburst behavior of RS Oph, and the recurrent 
nova phenomenon in general, has bcen reccntly discussed in [5]. 

The outbursts of PU Vul and other SyNe develop very slowly: the rise to 
máximum takes months to years, and their decline to the pre-outburst state can 
take decades. The outburst of PU Vul began in late 1977, and while the optical 
máximum was reached in 1979, the system has not yet recovered to its pre-outburst 
stage. The record-holder among the SyNe is however AG Peg: the constant- 
luminosity phase lasted for about 120 years, and the plateau luminosity, ~ 3500 ± 
1000 Lq, was the lowest among the whole sample (see Fig.6 in [T5]). Most of 
the SyNe develop an A/F-type spectrum at optical máximum, similar to that 
of classical novae; however without large violet displacemen ts o f the absorption 
lines which may indicate lack of massive optically thick wind [12]. In V1016 Cyg, 
V1329 Cyg, and HM Sge such a spectrum was not detected, possibly because of 
a lack of early cnough spectroscopic observations. 

The outburst evolution of PU Vul has been recently studied by Kato et al. 
[T2l [T3"j . Their analysis of the observed optical and UV light curves with quasi- 
evolution model consisting of a series of static solution resulted in the best fit 
mass of the white dwarf of ~ 0.6 M Q (FigQT])- Because of the low white dwarf 
mass, the TNR event is very slow and quiet without m assi ve optically thick wind, 
the white dwarf can retain most of the accreted mass [T^] . 

FigEDshows devclopment of the effective temperature and radius of the hot com- 
ponent of PU Vul, and its evolution in the HR diagram together with different 
models from ¡13] . The hot component evolves towards higher effective tempera- 
tures, and its radius shrinks as the nova declines from optical máximum. Thcsc 
changes are substantial; the effective temperature was ~ 160 000 K in 2003, while 
T e ñ ~ 6000-7000 K at visual máximum. Similar ly, the radius decrcased from its 
máximum valué of ~ 100 R© in 1982-83 to ~ 0.13 R© in 2003. Dcspitc thcsc 
temperature and radius changes, the luminosity remained constant at least till 
mid 2000's, in agreement with the timescales predicted by the best fit model with 
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. 5. Top: The evolution of the effective temperature and radius of PU Vul. Lines 
denote the blackbody photospheres with different mass-loss rates: 4 x 10 -7 (dash- 
dotted), 5 x 1CT 7 (thin solid), 6 x 1CT 7 (dotted), 7 x 1CT 7 (thick solid), and 8 x 1CT 7 
(dashed), respectively. Bottom: The HR diagram. Lines denote HR tracks for 
different models. Solid lines: 0.6 M Q , X=0.7, Z=0.01 (thick black - the best fit 
model), 0.6 M Q , X=0.5, Z=0.006 (magenta), 0.7 M Q , X=0.5, Z=0.006 (green), re- 
spectively. Dashed green line denotes a 0.7 M© with X=0.7, and Z=0.006. Dotted 
lines represent cold 0.5 Mq white dwarfs with different X and Z. The lowest dotted 
line represents a hot 0.5 Mq white dwarf with X=0.7, and Z=0.02. Arrows indicate 
evolution timescales for the best fit model. 
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0.6 M white dwarf [13]. 

The outburst behavior of PU Vul is similar to that of other SyNe ¡20l [T5l l23l [2Tj . 
although their plateau luminositics and timescales of the outburst significantly 
diffcr. Fig|6¡]presents correlation between thc lifctime (í p i) and luminosity (£ p i) 
of SyNe during their plateau phase. Apparently, the lower is the luminosity, the 
longer time the nova spcnds on the plateau. Such a behavior is expected from the 
thcoretical studies, and reflects dependence of both parameters on the white dwarf 
mass. For two of thesc novae, AG Peg and V1329 Cyg, the white dwarf masses 
have been derived from radial vclocity curves ¡16] , and they are in good agreemcnt 
with the theoretical prediction. AG Peg with its very low, 0.46 + / — 0.10 M©, mass 
of the white dwarf is the slowest nova even recorded whereas one of the fastest 
of the SyNe, V1329 Cyg, contains much more massivc, 0.74+ / — 0.08 M , white 
dwarf. 
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Fig. 6. The plateau lifetime- 
luminosity relationship for 
symbiotic novae. To as- 
sociate mass and luminos- 
ity on the top axis of 
the figure, the mass-lumi- 
nosity relationship for the 
plateau phase from 
was adopted. 



Interestingly, the presence or absence of the A/F-type spectrum at máximum 
is not obviously related to the plateau characteristics which may indicate that 
this is an observational effect rather than a real physical difference. In particular, 
V1329 Cyg, V1016 Cyg and HM Sge are the most luminous SyNe, and if their 
A-typc phase was shorter than in the other SyNe we could have lost it. Anothcr 
pos sibility is that they underwent non-dcgcncratc shell flashes on hot white d warf s 
[Ti] . Unfortunately, there are only few pre-outburst observations for the SyNe [2D] . 
On oíd objective prism plates, three of them, V1016 Cyg, V1329 Cyg, and PU Vul, 
look as M giants, and only V1016 Cyg shows a strong Ha emission indicating thc 
presence of a hot whit e dw arf prior to its eruption. Thc situation is much better 
in the case of RX Pup [2T| which have exhibited irregular photometric variability 
accompanied by significant spectroscopic changes. In particular, in 1941 is showed 
optical continuum with very strong high-excitation emission linos, whereas its first 
spectrum taken in 1894, resembled that of r¡ Car. The hot component clearly 
showed more activity before its 1970 nova outburst than any other symbiotic 
nova, and its white dwarf should have been the hottest among them. However, 
its evolution in thc HR diagram during its nova outburst is si mila r to that of PU 
Vul, with thc A/F-type spectrum maintained during 1972-79 ¡21) . 
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5. The Case of V407 Cyg 

The most recent nova eruption of V407 Cyg is both unique and challenging. 
The outburst started in the first week of March 2010, and its optical spectropho- 
tometric development is very similar to that of RS Oph. In particular, the very 
fast rise to máximum and decline within a few months to its pre-outburst visual 
brightness indicates that the white dwarf is massive as in the SyRNe (FigJT]). It 
is very surprising because V407 Cyg was hitherto known as a classical D-type 
symbiotic system which means that it must be wind-accreting system. 
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Fig. 7. The AAVSO optical light curve of V407 Cyg compared to that of the 2006 
outburst of RS Oph (insert). The horizontal scale (JD) is the same for both light 
curves. 

The Mira component of V407 Cyg is extreme, in the sense that it has the longest 
pulsation period among symbiotic Miras (P pu i = 763 days [27j). and the only one 
which is Li-rich [35] ■ Such a Li enhancement is very rare among evolved giants 
- it has been found in the LMC Miras with very long, > 400 days, periods, and 
only in one galactic Mira with a period of ~ 700 days - and it is usually explained 
as a consequence of hot bottom burning (HBB) ¡32] ■ Since the HBB oceurs only 
in stars with initial masses in the range of 4-8 Mg, the first formed white dwarf 
component of V407 Cyg should have had a massive progenitor, > 4-8M©, so it 
can be indeed very massive. 

There are also other remarkable similarities between V407 Cyg and the SyRNe. 
Bcforc the recent TNR outburst, V407 Cyg showed activity (high and low stages) 
similar to that observed in CH Cyg and other accretion-powered symbiotic stars 
¡2"T] , as well as the intrinsic variability of the hot components of the SyRNe be- 
tween their TNR nova outbursts. In particular, in all these systems the activity 
is characterized by very similar timescalcs, and a variable A/F-type shell source 
with ¿uv/opt ~ 10-1000 Lq accompanied by flickering, similar to that observed 
in cataclysmic variables, appears at the bright (high) stage Q] IS1 ES HH1 l2"Tj . 
Moreover, the same timescales characterize múltiple outburst activity of Z And- 
type symbiotic stars, and the same correlation between flickering and activity was 
found in Z And itself during its last series of eruptions [55]. The activity of Z 
And-type symbiotic stars and the high and low states of the SyRNe, CH Cyg, and 
other similar systems have been associated with unstable disc-accretion onto the 
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white dwarf (e.g. [15]). In the case of V407 Cyg, the hot component luminosity 

during the active state in 1998, ~ 5OOL [37], requires M acc ~ lO _7 M /yr (the 
distancc of ~ 2.7 kpc with the period-luminosity relation was adopted), similar to 
that expected in the SyRNe. Although such an accretion rate is somewhat high 
for a wind-accreting system, it is not impossible. Symbiotic Miras are known as 
sources of strong stcllar winds with rates of about of 10 -6 -10 -5 M /yr, and those 
with the longest pulsation periods have the highest mass-loss rates [5]. Moreover, 
recent theoretical simulations of interactions and wind-accretion in symbiotic bi- 
narles with dusty Mira donors show that once the dust-formation radius compares 
within a factor of a few to the Rochc-lobe radius, one can expect wind Rochc-lobe 
ovcrñow, whcre most of th e wi nd is trapped within the binary and can cfficiently 
accreted by the companion [22] . It is also very likely that the binary orbit is eccen- 
tric, and enhanced mass transfer rate near periastron passage causes a brightening 
of the accreting component. 

Finally, is V407 Cyg the first recurrent nova in D-type symbiotic system? The 
answer is NO because by definition at least two TNR cxplosions must be recorded 
whereas the 2010 outburst of V407 Cyg is the first one of this type. There have 
been some suggestions, also during this meeting, that the 1936 brightening was 
due to another nova outburst, however both amplitude and timescale pointto a 
high state similar to that in 1998, although somewhat fainter and shorter ¡27] . 

6. Summary 

Most of the symbiotic hot components are low mass white dwarfs powered by 
more or lcss stablc thermonuclear burning of the H-rich material accreted from 
their giant companion. However, they seem to be systcmatically cooler than the 
steady models predict. 

A thermonuclear nova explosión is a very rare phenomenon in symbiotic stars. 
Thcre are two kinds of symbiotic novae. The ordinary SyNc (9 hitherto recorded) 
are very slow and quiet (with no massive optically thick wind) TNR cxplosions 
on wind-accreting lowto medium-mass white dwarfs which can retain most of the 
accreted mass (e.g. [12j[l3]). The SyRNe are TNR explosions on very massive, 
^> 1.2M Q , white dwarfs accreting at very high rates, ~ 10~ 7 ' M Q /yr via Roche- 
Lobe overfiow. Low mass of the nova ejecta, ~ 1O _7 M , in RS Oph indicates 
that the white dwarf mass can grow, and thus they can be progenitors of the SN 
la. However, the most recent symbiotic nova, V407 Cyg, does fit well this picture 
because it seems to harbor a massive, >^ 1M , white dwarf as do the SyRNe. 
The white dwarf is accreting at relatively high rate, at least occasionally rcaching 
the highest rates in the SyRNe. 

Finally, prior to its nova explosión, V407 Cyg has been recognized as accretion- 
powercd symbiotic system sharing pronounced similarity with CH Cyg and R Aqr 
[2Tj . and even Mira itself [25]. They are several other symbiotic systems with 
less luminous hot components which can be accretion-powered. It is very likely, 
that they all have experienced, and will continué to experience symbiotic nova 
cxplosions. 
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